Biophysics and Physicobiology

Amylose contributes to starch crystallinity, but the stability of packed amylose double helices in water at elevated temperature remains insufficiently characterized. Here, we used molecular dynamics simulations to test whether chain length affects the short-timescale stability of A-type amylose oligomers in water. Six systems differing in chain length (6, 12, or 24 glucose units per chain) and oligomer size (isolated double strand or dodecamer of six double strands) were simulated, and five independent 1-s production runs were analyzed for each simulated condition. Oligomers with six glucose units showed structural collapse accompanied by increased water penetration. By contrast, dodecamers with 12 or 24 glucose units largely retained packed double-helical organization over the simulated timescale, although fraying was observed at their ends. These results indicate that chain length and lateral packing strongly affect the early structural response of amylose-like crystalline segments in hot water. The present simulations do not establish the ultimate fate of longer oligomers at longer timescales, but they identify a relative stability difference that is relevant to molecular interpretations of hydration-driven disordering in starch.

The high energy-efficiency of life is hard to understand only with classical physics. Many efforts have been made to study its mechanism based on quantum mechanics; the progress is nevertheless slow due to lack of experimental evidence with living cells. Here, combining experiments on cells, tissues and mitochondria with a theoretical model, we demonstrate a quantum state of mitochondria, which can be employed to modulate ATP production in living cells. We found an anomalous 71.0-THz oscillation mode only in living cells and tissues, which is highly determined by intact structure of mitochondria, and cannot be assigned to any specific molecules. Based on experimental data, a quantum model of light-matter coupling was introduced to trace the origin of this mode. Our calculations suggest a quantum superposition state of functional mitochondrion that forms by the coupling of light and lipid CH2 bonds in functional cristae, and induces a splitting of the intrinsic CH2 vibration mode of 87 THz to two levels at 71 THz and 103 THz, respectively. The former can be observed only in living cells and tissues; whereas the latter falls in the range (90-110 THz) of biomolecular and water vibrations, thus indistinguishable. Additional experiments revealed this mitochondrial quantum state able to serve as an efficient channel to modulate ATP production. Our findings provide a quantum mechanics view for understanding living cells, and it will be interesting to further explore whether such quantum state could act as a channel for energy metabolism, and even information transmission in life.

G protein-coupled receptors (GPCRs) represent the largest membrane protein family in the human body and the major targets of FDA-approved medications. Unlike some other physiological proteins, GPCRs are highly plastic and often exist as collections of conformational ensembles in the environment to exert their function. Therefore, understanding how a GPCR adopts different conformational ensembles in response to ligands and signaling partners is essential for drug development. 19F-NMR has emerged as one of the most promising tools in this regard; however, a well-established 19F-labeling system is still lacking, including an effective 19F probe that meets both environmental susceptibility and NMR sensitivity requirements, even though progress has been made dramatically in recent years. Thus, this manuscript reports a novel conformational probe (4-(bromomethyl)-ethoxy-2-(trifluoromethyl)benzene, BTMB) that offers paramount benefits for future GPCR conformational dynamics studies. To rationalize the mechanism behind BTMB's superiority for future probe development, we performed MD simulations, referenced to the current state-of-the-art probe, BTFMA. Our study suggests that the spatial proximity of aromatic residues significantly influences the chemical shift dispersion and linewidths of related 19F-probe-profiled conformational states, as well as the rigidity and compactness of -CF3 in the target protein.

G protein-coupled receptors (GPCRs) are important targets for drug discovery because they are the largest and most diverse family of membrane proteins in the human body. As exemplified by {beta}2-adrenergic receptor ({beta}2AR), which is a typical class-A GPCRs, they are prone to denaturation and inactivation after solubilization. Although diverse techniques have been developed, current structural and functional analyses are still limited to proteins with relatively stable and high expression. However, some mutant variants and misfolded proteins involved in diseases have extremely low expression levels and may be difficult to analyze. To overcome these limitations, we established a novel analytical platform for evaluating the ligand-binding ability of full-length {beta}2AR as a model at the single-molecule level without purification and addressing challenges such as membrane proteins with low expression levels and structural instability. This method enables the direct use of unpurified receptors immediately after solubilization and allows us to use only a small amount of sample ([~]10 ng) for observation and to distinguish between specific and nonspecific ligand binding by fitting. Furthermore, we unveiled the physical properties of detergents on the structural stability of solubilized receptors and found that the lateral pressure within detergent micelles affects ligand-binding ability. Detergents that provided a fluid microenvironment were able to maintain ligand-binding ability for several days even after solubilization; conversely, detergents that provided a rigid microenvironment caused the protein to lose its activity earlier. Our method could be a promising tool for the structural and functional analysis of membrane proteins untargeted until now. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/740883v1_ufig1.gif" ALT="Figure 1"> View larger version (46K): org.highwire.dtl.DTLVardef@1ca0ed1org.highwire.dtl.DTLVardef@92e757org.highwire.dtl.DTLVardef@9223edorg.highwire.dtl.DTLVardef@74a100_HPS_FORMAT_FIGEXP M_FIG C_FIG

Interleukin-13 (IL-13) is an immunomodulatory cell signaling cytokine that has been implicated in neurodegenerative disease and chronic inflammation. IL-13 binds to its low and high affinity receptors, IL-13 receptor 1 (IL-13R1) and IL-13 receptor 2 (IL-13R2), respectively, with residence times that vary accordingly. As the binding kinetics of the cytokine-receptor complexes influence cellular responses, we employed the molecular dynamics (MD) simulation-based{tau} -random acceleration molecular dynamics method ({tau}RAMD) to compute relative residence times for wild-type (WT) IL-13 and 19 IL-13 mutants to the two receptors. Comparison with experimental kinetic data shows that the{tau} RAMD computations capture the trends in residence times. Analysis of simulated dissociation trajectories of the cytokine-receptor complexes reveals two distinct dissociation pathways of IL-13 from each of the receptors. This study thus pinpoints key determinants of the interaction of IL-13 with its receptors which could be targeted for therapeutic design. Statement of SignificanceCytokines are regulatory proteins that bind to cell surface receptors and thereby send signals to the cellular interior. Interleukin-13 (IL-13) is a cytokine that has a low and a high affinity receptor. It has important physiological roles, and its deregulation is involved in diseases such as atopic dermatitis and asthma. We employed a molecular dynamics simulation-based method to compute the effects of changes in the sequence of IL-13 on the lifetimes of complexes of IL-13 and its receptors. Comparison with experiments supports the validity of the computational approach and analysis of the simulations reveals two distinct ways in which IL-13 dissociates from each receptor. These results thus provide a map for targeting IL-13 - receptor interactions for the design of therapeutics.

Imaging spectroscopy often requires expensive and complex equipment. Here we show a simple procedure for attaching a standard miniature fiber spectrometer to a conventional microscope, allowing easy integration of 2D imaging with 1D high-resolution spectral measurements. This combination provides much of the benefit of a full imaging spectrometer without the large equipment investment, and we provide instructions for modifying microscopes to this setup and the present measurements of living cells that demonstrate their performance. Using this setup, we compare the quantitative measurement of chlorophyll concentration in Physcomitrium patens moss using color imaging and spectral sampling.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

The mechanical properties of red blood cell (RBC) membranes are critical to their function in oxygen delivery, and changes to these properties as RBCs age affect their journey around the circulatory system, including clearance by the spleen. Such changes can also have significant health effects,including in cardiovascular disease and on the interactions between RBCs and malaria parasites. The function of RBCs requires them to be very soft, which together with the small size of the cells brings the energy required for measurable deformation of the membrane into the range of typical thermal energies. Red Blood Cells can be observed to flicker under optical microscopy, and these shape fluctuations can be quantified and used to obtain key biophysical parameters such as the tension and bending modulus of the membrane. Typically the shape of the cell's equator is extracted, and the mean power spectrum is obtained by time-averaging the power present in the normal modes of the thermal fluctuations. Flickering spectroscopy has been used extensively on RBCs, but has so far been very low throughput and with some technical limitations. Here we address issues related to active versus passive fluctuations, focusing and optics, camera exposure and sampling, and fast contour detection. In combination with an automated imaging system it is possible to measure thousands of cells in one day, with no user input. We validate this new pipeline by chemical modification of RBC mechanics, and by comparison with simulated fluctuations including each confounding effect. The methods and codes for this robust flickering analysis will allow consistent measurements across labs.

The development of NIR fluorescent dyes with favorable optical properties is of great importance especially for applications within the field of optical imaging. In this study, we report the design, synthesis and photophysical characterization of novel BODIPY-Xanthene hybrid dyes obtained via Knoevenagel condensation between an asymmetrical BODIPY core and the functionalized xanthene aldehydes intermediates. The synthesized dyes ex-hibit absorption (667 nm) and emission bands (720-722 nm) within the NIR region. Detailed UV-Vis and fluo-rescence studies were performed and molar extinction coefficients and stokes shifts also calculated. Moreover, findings highlighting the influence of structural modifications on the optical behavior of the hybrid dyes and their potential utility in optical imaging applications are reported.

Boron neutron capture therapy (BNCT) is one of the promising treatment methods for cancer. BNCT is based on the unique high ability of the non-radioactive boron-10 (10B) nucleus to absorb thermal neutrons. The absorption of a neutron by boron results in a nuclear reaction 10B(n,)7Li with 84% of the energy being released within the cell, thereby inducing cancer cell death. We propose a novel approach based on using lithium instead of boron in neutron capture therapy: the neutron capture reaction cross-section for lithium is 4 times smaller than that of boron, while the energy release is 2 times higher, and the primary advantage of lithium is that 100% of the energy is released inside the cell. A series of in vivo neutron irradiation experiments using the murine B16 skin melanoma model demonstrated the efficacy of lithium neutron capture therapy (LiNCT). Higher lithium concentrations in the tumor achieved through intraperitoneal administration, compared to the oral route, contributed to a significant increase in animal survival and a significant decrease in tumor growth. The results of our study provide entirely new opportunities for the development of neutron capture therapy, both with lithium and with a possible lithium-boron combination.

It is important to determine whether native state dynamics drive the misfolding and oligomerization of the prion protein, which are important events in prion disease, and how they are modulated by conformational conversion. Native (N) mouse prion protein (moPrP) is known to form small (OS) and large (OL) oligomers rich in {beta}-sheet, and in this study, photoinduced electron transfer-fluorescence correlation spectroscopy (PET-FCS) has been used to characterize intramolecular dynamics within individual monomeric units in both isolated OS and OL, as well as the diffusion properties of the oligomers. It is estimated that OS and OL comprise of about 15 and 55 monomeric units, respectively. Microsecond dynamics at each of the two regions that are the 1-3 and 2-3 interfaces of native protein are distinct in N, OS and OL, although they occur on very similar timescales. Analysis of the evolution of the distribution of diffusion times, determined using the maximum entropy method, indicates heterogeneity in the oligomerization reaction. Analysis of the change in the fluctuations which occur in two different timescales in the native state ensemble shows that they are damped more at the erstwhile 1-3 interface than the erstwhile 2-3 interface. The difference in the extent of damping at the erstwhile 1-3 and 2-3 interfaces can be explained on the basis of the structural changes known to occur across each region. The changes in dynamics occur concurrently in both regions, indicating that the structural changes accompanying conformational conversion also occur simultaneously during the oligomerization of moPrP.

The process of drug unbinding is of immense importance in the field of biophysics and therapeutics. The behavior of these systems is greatly influenced by their thermodynamic and kinetic properties. Therefore, it is crucial to accurately estimate the ligand binding free energies and rate of ligand dissociation, yet these processes are often governed by rare event transitions that lie beyond the reach of standard brute-force molecular dynamics simulations. While enhanced sampling simulations offer a solution, their efficacy is strictly contingent upon the selection of appropriate collective variables (CVs) which is non-trivial for complex systems like protein-ligand complexes. In this study, we present a method to derive optimized CV from transition state region (TS) via an interpretable machine learning (ML) model, Elastic Net. By employing some physically intuitive order parameters, the derived optimized CV from the TS-region greatly accelerate ligand binding-unbinding transitions and achieves rapid free energy surface (FES) convergence across diverse systems including buried and solvent exposed active sites such as Trpsin-benzamidine complex, host-guest systems and sodium epoxidase etc. Intriguingly significant contribution of the ligand hydration is found in the optimized CV which depicts crucial role of solvent in driving ligand binding-unbinding transitions. The estimated binding free energies for different protein-ligand complexes match quite well with experiments, while maintaining a low computational cost. The derived optimized CV is also used to calculate the ligand residence times across different systems and calculated residence times are within the experimental range for all systems, again with very little computational costs. Moreover, we show that the optimized CV constructed from TS region via an interpretable ML model is transferable across diverse systems, offering a robust and scalable framework for drug discovery and investigation of complex biomolecular recognition.

Far-Red (FR) Light Photoacclimation (FaRLiP) enables cyanobacteria to extend photosynthetic activity into the far-red region by extensively remodelling Photosystem I (PSI), including the replacement of several core subunits with paralogs that coordinate the red-shifted chlorophyll f (Chl f). The binding positions of Chls f are still a matter of debate, with the most recent structural findings supporting the location of a single Chl f molecule within the reaction centre (RC) at the so-called A-1B site. This was in turn suggested to strongly affect electron transfer (ET) directionality leading to an almost monodirectional transfer along the B branch in FR-PSI RC. Here, we directly probe ET in FR-PSI by characterising the photogenerated [P700A1-] spin-correlated radical pair using complementary pulse and Time-Resolved (TR) Electron Paramagnetic Resonance (EPR) spectroscopy at cryogenic temperature. Electron spin-echo decay kinetics are distinctly biexponential, indicating the formation of two charge-separated states. Consistently, out-of-phase ESEEM traces are quantitatively described by two modulation frequencies arising from different dipolar interactions, while TR-EPR spectra are accurately simulated by the combined contributions of [P700A1A-] and [P700A1B-] radical pairs. These results provide direct spectroscopic evidence that both the A and B branches remain photochemically active in FR-PSI. The conservation of bidirectional ET, even when considering the presence of a single Chl f molecule in the RC, further implies that the two radical pairs originate from a common primary electron donor. This finding identifies P700 as the most likely primary donor and argues against a mechanism in which the RC Chl f initiates charge separation.

A universal mechanism in molecular evolution is functional and structural divergence of members of a protein family. The ability of AlphaFold to predict atomic-resolution protein structures promises to accelerate insights into this process. We study the interplay of changes in sequence, structure, and function in photoactive yellow protein (PYP), a family of bacterial blue light photoreceptors. Halorhodospira halophila contains two PYP homologs that diverged to 60% sequence identity, differ 100-fold in the lifetime ({tau}pB) of their pB signaling intermediate, and display altered peak wavelengths ({lambda}max) for color sensing. We resurrected ancestral PYPs and determined these properties along the resulting recapitulating evolutionary divergence. The resurrected ancestral PYP is functionally similar to PYP1, indicating divergence on the path to PYP2. AlphaFold predictions for PYP2 and these ancestral proteins revealed the absence of structural changes compared to the crystal structure of PYP1. To experimentally validate these predictions, we optimized second-derivative Fourier transform infrared (FTIR) spectroscopy. The FTIR spectra of PYP1 and 2 and their resurrected ancestral proteins demonstrated clear differences in their secondary structure. These results demonstrate an important limitation of AlphaFold and show how ancestral sequence reconstruction combined with spectroscopic approaches yields insights into divergence in a protein family.

Phosphatase and tensin homolog (PTEN)-induced putative kinase 1 (PINK1), a key regulator of mitophagy, has been linked to the pathogenesis of Parkinson's disease (PD). PINK1 recruits Parkin, an E3 ubiquitin ligase, triggering mitophagy in response to mitochondrial damage. During mitophagy, the quantity, stability, and activity of PINK1 must be strictly regulated; however, the mechanisms governing these parameters under cellular stress are still unclear. Herein, we determined the structural basis for PINK1 maturation mediated by heat shock protein 90/cell division cycle 37/FK506-binding protein 51 (HSP90/CDC37/FKBP51) chaperone complex. We identified PINK1-associated proteins using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and determined the structures of the complexes using Cryo-Electron Microscopy (Cryo-EM). Results showed that FKBP51 potentially interacts with a conserved leucine-proline-phenylalanine (LPF) motif on the activation loop of PINK1 and negatively regulates PINK1 functions in mitophagy. A PINK1 mutation located at the FKBP51 recognition site is linked to mitophagy deficiency, which can be partially rescued by specific inhibition of FKBP51. These findings reveal a general mechanism for PINK1 recognition by the HSP90/CDC37/FKBP51 chaperone complex and suggest a potential approach for upregulating PINK1 activity, which is impaired in PD.

Genome scale metabolic models provide detailed mechanistic descriptions of cellular metabolism, whereas black box growth models capture physiological behaviors using a small number of effective parameters. However, the quantitative relationship between these two modeling scales remains unclear. In this study, we develop a thermodynamic framework that connects black box growth models with thermodynamically constrained genome scale metabolic models. By linking black box model parameters with Gibbs energy dissipation rates derived from genome scale metabolic models, we demonstrate that coarse grained physiological descriptions can be obtained from detailed metabolic networks while preserving their thermodynamic foundation. We validate this framework in both Escherichia coli and yeast, showing that the resulting black box models reproduce key physiological behaviors, including growth dynamics, biomass yield, and overflow metabolism observed experimentally. Our results indicate that black box growth models and genome scale metabolic models are connected through shared thermodynamic constraints, revealing a consistent thermodynamic basis across different levels of metabolic description. This framework provides a general approach for integrating detailed metabolic networks with simple black box growth models and enables efficient multiscale modeling of cellular metabolism.

Cholesterol is a key component of eukaryotic cell membranes, promoting membrane stability and modulating the function of many membrane proteins, including ion channels. In our previous work using purified human voltage-gated proton channel proteins, we showed that cholesterol inhibits the hHv1 channel by altering the conformational dynamics of its S4 segment, the key element that senses membrane voltage to control proton permeation. In the present work, we examined the effects of cholesterol analogs and potential sites in the hHv1 channel mediating cholesterol inhibition using site-directed mutagenesis and docking simulations. Our results showed that desmosterol, the immediate precursor of cholesterol, markedly attenuates cholesterol inhibition. Using single-molecule Fluorescence Resonance Energy Transfer (smFRET), we showed that desmosterol attenuates cholesterol inhibition by promoting the intermediate and open state conformations of the S4 segment. Moreover, we identified multiple residues in the hHv1 channel that are critical for cholesterol inhibition, including Y141A in the S2 segment, which reduces cholesterol inhibition by nearly 3-fold. Our smFRET results showed that the Y141A mutation promotes the intermediate conformation in the S4 segment, which underlies the attenuation of cholesterol inhibition. Consistently, docking simulations also revealed multiple residues spanning the transmembrane domain, rather than clustered within a single localized pocket. Our work identified the key molecular determinant in the hHv1 channel that mediates cholesterol inhibition and also provided a mechanism linking the conversion between demosterol and cholesterol by DHCR24 to pH homeostasis in many cells, such as phagocytes, cardiomyocytes, neurons and microglial cells.

Imipramine is known as an older generation tricyclic antidepressant drug. It has been identified in prior studies that imipramine blocks Apolipoprotein E4 (ApoE4)-induced amyloid-{beta}(A{beta}) aggregation and is associated with an improved AD diagnosis [Johnson et al. Alzheimers Research Therapy, 2022, 14, 88]. Using NMR methods such as 1H-1H NOESY and Saturation Transfer Difference Spectroscopy, we demonstrate the binding of A{beta} monomers to imipramine when the full-length A{beta} (1-42) sequence is considered. The more abundant but less toxic form, A{beta} (1-40) does not show interaction with imipramine.

Amyloid {beta}-protein (A{beta}), which is known for its role in Alzheimer's disease (AD), is also involved in several physiological functions. Of the two predominant isoforms, A{beta}40 and A{beta}42, the latter is more aggregation-prone and more active in health and disease. In the brain, A{beta} interacts with free lipids, which may affect its assembly pathways. Here, thioflavin T (ThT) fluorescence spectroscopy is combined with molecular dynamics (MD) to investigate A{beta}40 and A{beta}42 co-assembly with free dimyristoylphosphatidylcholine (DMPC) and DMPC/cholesterol (CHOL) lipids. ThT measurements at several protein-to-lipid (PL) ratios show isoform- and lipid composition-dependent kinetics of A{beta} fibril formation with clear isoform-dependent lipid effects at a PL ratio of 1:5. Early co-assembly stages of six A{beta}40/A{beta}42 with DMPC and DMPC/CHOL lipids at this PL ratio are then examined by MD. MD results indicate that lipids alter the secondary structure more and interact more avidly with A{beta}40 than with A{beta}42. In turn, A{beta}40 increases the planar order of lipid bicelles more than A{beta}42. Unlike A{beta}40, A{beta}42 decreases the nematic order of CHOL. In DMPC/CHOL lipids, A{beta}42 but not A{beta}40 exhibits two A{beta}-lipid interaction modes associated with mutually exclusive A{beta}42-lipid and A{beta}42-A{beta}42 interactions. Our MD findings provide a plausible explanation for an experimentally observed increase in the lag time for A{beta}40 fibril formation in DMPC lipids and even more so in DMPC/CHOL lipids that contrasts with much smaller effects of lipids on A{beta}42 fibril formation and offer additional atomistic-level insights into A{beta} co-assembly with lipids amenable to experimental testing.

Eukaryotic translation initiation is tightly regulated by interactions among translation initiation factors (eIFs) that ensure accurate start codon selection. The translation regulator, eIF5 mimic protein 1 (5MP1) contributes to this process by competing with eIF5 for binding to eIF2, thereby increasing the stringency of translation initiation. Despite its important regulatory role and emerging involvement in tumorigenesis, structural information on human 5MP1 remains limited. Here, we report the near-complete backbone and partial side-chain NMR resonance assignments of the C-terminal domain of human 5MP1 (residues 250-419), carrying a W404E substitution that disrupts dimerization. The WT protein forms a dimer at NMR concentrations, which increases the effective size of the protein and also causes disappearance of peaks corresponding to aminoacids at the dimer interface due to conformational exchange. Backbone resonance assignments were completed for 96.4% of the non-proline residues. Secondary structure was analyzed using Chemical Shift Index (CSI) and compared with the AlphaFold structural model. Regions of disagreement between the experimental and computational secondary structure assignments were further examined using 15N-NOESY-HSQC spectra, allowing experimental validation of local structural features. While the AlphaFold model accurately reproduces the overall fold of the 5MP1 C-terminal domain, several localized discrepancies were identified, particularly near the N- and C-terminal regions of the domain, where experimental NMR data support alternative secondary structure assignments. These resonance assignments and experimentally validated structural features provide a foundation for future investigations of the molecular interactions, dynamics, and functions of 5MP1 in translation initiation.